Correlations in Quantum Physics 



Ross DorneiQ 

Blackett Laboratory, Department of Physics, Imperial College London, London SW7 2AZ, U.K. and 
Clarendon Laboratory, Department of Physics, University of Oxford, Parks Road, Oxford 0X1 3PU, U.K. 

Vlatko Vedral 

Clarendon Laboratory, Department of Physics, University of Oxford, Parks Road, Oxford 0X1 3PU, U.K., 
Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, 117543, Singapore and 
Department of Physics, National University of Singapore, 2 Science Drive 3, 117542, Singapore 

We provide an historical perspective of how the notion of correlations has evolved within quantum 
physics. We begin by reviewing Shannon's information theory and its first application in quantum 
physics, due to Everett, in explaining the information conveyed during a quantum measurement. 
This naturally leads us to Lindblad's information theoretic analysis of quantum measurements and 
his emphasis of the difference between the classical and quantum mutual information. After briefly 
summarising the quantification of entanglement using these and related ideas, we arrive at the con- 
cept of quantum discord that naturally captures the boundary between entanglement and classical 
correlations. Finally we discuss possible links between discord and the generation of correlations in 
<•—{ thermodynamic transformations of coupled harmonic oscillators. 
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I. INTRODUCTION 

Correlations play a prominent role in many-body physics, statistical physics and information technology. In this 
article, however, we focus on another motivation for studying correlations, namely that they offer the best way of 
understanding the key differences differences between quantum and classical physics. Quantum mechanics allows 
for the full range of correlations between events. There is a maximal correlation between earlier and later states 
^ , of quantum systems undergoing a unitary evolution. Also, if a measurement is repeated successively on the same 
system, it will always give the same outcome. At the other extreme lie completely uncorrelated quantum events, 
such as measuring a spin half particle successively in two complementary basis. All other partially correlated events 
between maximal and zero correlation can also be found in quantum physics. Firstly, we will explain how information 
theory can be used to quantify correlations in general. We then show how this was applied by Everett and Lindblad to 
elucidate the quantum measurement process. This is followed by a brief summary of the two main types of quantum 
correlation, entanglement and quantum discord. Finally, the role of correlations is explored in the field of quantum 
thermodynamics, where we explore links between correlations and the amount of work done by two couple harmonic 
oscillators. 
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II. SHANNON'S INFORMATION THEORY 

> . 

• —i Given two random variables, X with set of possible outcomes {x{\ and Y with the set of possible outcomes {yj}, 

we define the joint probability distribution P(X = Xi,Y = yj) = p(xi,yj) with corresponding marginal probability 
distributions P(A = at) = p{xi) = J^ - p{xi,yj) and P(Y = yj) = p(yj) = ^2iP(xi,yj). The degree of correlation 
between the random variables X and Y is encoded within their joint probability distribution and best quantified 
within the framework of Shannon's information theory [lj. To that end, we begin our discussion by introducing the 
Shannon entropy of the random variable X, defined as 

S(X) = S(p(x)) = -^pixjlnpixi), (1) 

i 

and analogously S(Y) — —^2jP{yj)^p(yj). As is well known, the Shannon entropy provides a measure of the 
uncertainty in a given random variable. By straightforward generalisation of Eq. ([!]) the uncertainty in the joint 
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distribution of X and Y is defined by 

S(X,Y) = S(p(x,y)) = -^p{x l ,y, J )\np{x l ,y j ). 

From this, we are able define the the Shannon mutual information, 

I S (X : Y) = S(p(x)) + S(p(y)) - S(p(x, y)), (2) 

This quantity, being positive, provides us with an intuitive interpretation of correlation between X and Y. Namely, 
that the sum of the uncertainties in p(x) and p(y) alone is greater than the uncertainty in the their joint distribution 
p(x,y). Next, by introducing the Shannon relative entropy, 



S(X\\Y) = S( P (x)\\p(y)) In 



p{y 3 y 



h3 

the Shannon mutual information Eq. ([2| can be rewritten in the equivalent form 

I s (A:B) = S(p(x,y)\\p(x)xp(y)). 

In this sense, the mutual information represents a distance between the joint probability distribution p(x,y) and the 
product of its marginals p{x) x p(y). As such, it is intuitively clear that Is is a good measure of correlations, since it 
shows how far a joint distribution is from the corresponding product distribution in which all the correlations have 
been destroyed. Let us now view this from another perspective. Suppose that we wish to know the probability of 
observing the event yj given that the event Xi has been observed. The relevant quantity is the conditional probability, 
given by 

P{Y = y J \X = x l )^p x ^y 1 )- P{x ^ ) 



P(Xi) 

This motivates us to introduce the conditional entropy, Sx{Y), as 

Sx(Y) = -^2p(xi)^2p Xi (yj)lnp Xi (yj) 

i 3 
i,3 

The conditional entropy tells us how uncertain we are about the value of Y once we have learned the value of X. Now 
the Shannon mutual information Eq. ([2| can be rewritten once more as 

I S (X : Y) = S(Y) - S X (Y) = S(X) - S Y (X). (3) 

Thus, as its name indicates, the Shannon mutual information measures the quantity of information conveyed about 
the random variable X (Y) through measurements of the random variable Y (X). As this quantity is positive, it tells 
us that the initial uncertainty in Y (X) can in no way be increased by making observations on X (Y). Note also that, 
unlike the Shannon relative entropy, the Shannon mutual information is symmetric i.e., Is(X : Y) = Is(Y : X). 

Shannon's original motivation for developing this framework was the practical need of quantifying the channel ca- 
pacity of classical communication channels. In this scenario, X and Y represent the message sent down and received 
from the channel respectively. Given that process of performing a measurement can be viewed as a communication 
channel between the states of the system and the observer, it was only natural that Everett would use Shannon's infor- 
mation theory to characterise quantum measurements within his relative state interpretation of quantum physics [2]. 
The novelty of Everett's approach was that it became possible to discuss the information in quantum measurements 
without recourse the Born projection postulate. In this view, a quantum measurement is nothing but the establish- 
ment of correlations between the system and the measuring device. Everett did not discriminate between quantum 
and classical correlations, though he implicitly worked only with classical correlations as we now understand them. 



III. EVERETT'S VIEW OF CORRELATIONS 



According to the relative state view of quantum physics, a quantum measurement never gives conclusive results, 
but instead establishes correlations between the system and the measuring device. As an example, for a system of 
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a single qubit a typical Everettian style measurement therefore takes initially uncorrelated states of the system and 
measuring device to the post-measurement, correlated state 

|VA B )=«|0)|m 1 )+/3|l)|m 2 ), (4) 

where |0) and |1) are orthogonal states of the system, \m\) and \rri2} are states of the measuring device with overlap 
(60 1 61) = e and \a\ 2 + \(3\ 2 = 1. Intuitively, the degree of correlation between the system and measurement device 
depends upon e. This intuition can be quantified within the framework of Shannon's information theory by first 
introducing one observable pertaining to the system X = xi\xi){xi\ + x-i \x^ (a^l, and one observable pertaining 
to the measuring device Y = yi\yi)(yi\ + J/2 lz/2) (2/2 1 - The outcomes of measurements of these observables give rise 
to the two probability distributions, p{xi) = \(4>AB\xi)\ 2 an d p(yj) = | (V'Aslyj) | 2 as well as the joint distribution 
p(xi,yj) = I (■z/'ab l^i) I J/j ) 1 2 • With this we can straightforwardly quantify the uncertainty in a given observable by the 
Shannon entropy Eq. ([I]) and the correlation between the observables X and Y via the Shannon mutual information 
Eq. (2). 

At this point, we make the important remark that while Everett showed that Shannon's information theory is 
capable of quantifying some correlations in quantum theory, it does not tell us the whole story; While the Shannon 
mutual information can tell us about correlations between the observables X and Y , the correlation in post measure- 
ment quantum state itself Eq. Q, for instance, must be quantified using the quantum generalisation of Shannon's 
information theory. This will be the focus of Sees. TVffVI 



Everett went on to establish the properties that any must be possessed by any 'good' measure of correlations between 
two random variables. Namely, he proved that if either or both of the variables undergo a local stochastic evolution, 
then the amount of correlations between them can not increase (in fact it usually decreases, see also Ref. }3i). As far as 
quantum measurements are concerned, this describes the physically intuitive fact that correlations established during 
the quantum measurements cannot be increased, i.e., the measurement cannot be made to yield more information, by 
further processing the system and the apparatus separately. We note however, that when local stochastic processes 
are correlated they virtually become global, and therefore the correlations between the systems can increase as well 
as decrease. A generalisation of this logic is also central to the theory of measuring quantum correlations and the 
paradigm of local operations and classical communication (LOCC). We will return to this in Sec. [V] Next, however, 
we discuss another important step in the development of correlations in quantum physics, due to Lindblad. 



IV. LINDBLAD'S ANALYSIS OF QUANTUM MEASUREMENTS 



Lindblad built on Everett's initial work by generalising his theory of quantum measurements to encompass mixed 
states in contact with an environment !4 . In doing so, Lindblad employed the quantum versions of many of quantities 
from Shannon's information theory. The formulation of classical probability theory that is most naturally generalized 
to quantum states is provided by Kolmogorov [5], and an exposition expressing similarities with von Neumann's 
Hilbert Space formulation [H] due to Mackey [7] (see also Ref. [5]). For our purposes, however, it suffices to introduce 
the relevant quantities ad hoc by straightforward analogy with their classical counterparts. 

First, for a quantum system described by a density matrix p, we define its von Neumann entropy by 

Sn(p) = -Tr[plnp], (5) 

This can be considered as the proper quantum analogue of the Shannon entropy Eq. (JlJ [5]. For an observable A of 
the system described by p, the Shannon entropy S(A), describing the uncertainties in the values of the observables, 
is equal to the von Neumann entropy S'n(p) only when [^4, p] = i.e., A is a Schmidt observable. Otherwise, 

S(A) > S N (p). 

The physical interpretation of this result is that there is more uncertainty in a single observable than in the whole of 
the state, a fact which entirely contradicts our expectation. 

The concept of mutual information Eq. |2]) can also be easily extended to quantum systems |10j . For a bipartite state 
Pab and corresponding reduced states pa = Ttb [pab] and pb — Tta \pab\i t ne von Neumann mutual information 
between the two subsystems is defined as 

In(pa ■ Pb;pab) = S n (pa) + S N {p B ) - S n (pab), (6) 

which refers to the correlation between the whole subsystems rather than relating two observables only. Following 
the Shannon mutual information, this quantity can be interpreted as a distance between two quantum states by 
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introducing the von Neumann relative entropy, (in fact, this quantity was first considered by Umegaki [TT], but for 
consistency reasons we name it after von Neumann) . The von Neumann relative entropy between two quantum states 
a and p is defined as 

S N (a\\p) =Tr[<7(ln<7-lnp)]. (7) 

Now, the von Neumann mutual information Eq. ^ can be understood as a distance of the state pab to the uncorrelated 
state pa® Pb, 

In(pa ■ Pb;Pab) = S n (pab\\pa ® Pb)- (8) 

Lindblad's main contribution was to calculate and compare the classical and quantum mutual information during 
a measurement. This analysis was motivated by thermodynamical considerations. A question he considered is how 
much entropy is generated during a quantum measurement. The entropy increase is fundamentally linked with the 
irreversibility of the quantum measurement and the correlations created between the system and an apparatus. We 
will return to the connection between correlations and thermodynamics in Sec. |VII| 

We now discuss a relation concerning the entropies of two subsystems. One part of it is somewhat analogous to 
its classical counterpart, but instead of referring to observables it is related to the two states. The relation is the 
Araki-Lieb inequality |12) . 

S-n(pa) + Sn{pb) > S-m(p) > \Sn(pa) - Sn(pb)\, (9) 

which can be considered one of the most important results in the quantum theory of correlations. Physically, Eq. |9]) 
implies that their is more information (less uncertainty) in an entangled state than if the two states are treated 
separately. One the one hand, this arises naturally since by treating the subsystems separately we have neglected the 
correlations between them. However, to appreciate the extent to which this is a counter-intuitive result we consider 
the following example. Suppose a two level atom is interacting with a single mode of an EM field as in the Jaynes- 
Cummings model [13]. If the overall state is initially pure, and the whole system is isolated then the entropies of the 
atom and the field are equally uncertain at all the times. But this is not expected since the atom has only two degrees 
of freedom and the field infinitely many! This, is possible however, as by the second observation, two level atom, is 
only entangled with two dimensions of the field. 



V. QUANTIFYING ENTANGLEMENT 



The global state of a composite system is entangled if it may not be written as a product of the individual subsystems 
[T4l [T5] . This implies that though the global state of the composite system may be known with certainty, the state of 
each individual subsystem is uncertain, or, as Schrodinger put it, "The best possible knowledge of a whole does not 
include the best possible knowledge of its parts" (THj (c.f. Eq. (|9|). This reasoning can be quantified by considering 
the two qubit state 

\*ab) =coB-|0A0 B ) + Bin-|l i ilB>, (10) 

which, following the above definition is entangled for < 9 < n. For pure bipartite states, the amount of entanglement 
present is uniquely quantified by the von Neumann entropy Eq. (5]) of (either of) the subsystems [171 . Though the 
von Neumann entropy of the global state is zero, that of each of the reduced states is finite and determined by 9, 

9 9 9 9 

Sn(pa) = S n (pb) = - cos 2 - In cos 2 - - sin 2 -lnsin 2 -. (11) 

The entanglement is thus the lack of information concerning the reduced states of the subsystems in light of maximal 
information regarding the global state. Conversely, the more entangled the global state is the less information we 
have concerning the states of the subsystems. 

For a seperable state (9 = 0,7r), the von Neumann entropy of the global and reduced states are identically zero. 



Conversely, when cos(6>/2) = sin(0/2) = l/v2) the state Eq. (10 1 is maximally entangled and its reduced states are 
maximally mixed i.e., S'n(pa) = 'S'n(pb) — hi 2. It follows that the von Neumann mutual information Eq. ([6]) of a 
maximally entangled state /n(pa;/?b : Pab) = 2/n2, is twice the maximum possible value of the Shannon mutual 
information between two random variables. Lindblad was the first to attribute this stronger-than-classical correlation 
to entanglement, though this assertion is only accurate for pure bipartite states and in general, entanglement does 
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not encapsulate all forms of non-classical correlation. Note that for other values of < 9 < ir, the state is entangled 
but the von Neumann mutual information may not exceed the maximum possible Shanon mutual information. 

More generally, for finite dimensional N-partite states, entanglement presents any form of correlation that cannot 
be captured by a separable state of the form 

a = ^p{i)n ( ( ) ® ...®t:<$, (12) 

i 

where 7r„ is the reduced density matrix of the ./V th subsystem. Operationally speaking, seperable states are those 
that are preparable under the paradigm of LOCC, while all other states are entangled [HI [19]. The entanglement in 
a quantum state can thus be quantified by its distance to the closest separable state [TTJ, 120] • This distance can be 
expressed in a number of ways but the related details will not trouble us at present (see e.g. Ref. [21]). For infinite 
dimensional (continuous variable) systems, the theory of entanglement for the physically relevant set of Gaussian 
states is also well developed [221 [22] . 

Entanglement aside, there are other forms of correlation in quantum systems with no counterpart in classical 
information theory. Just as with entanglement, the degree of other quantum correlations in a given state may be 
measured by its closest distance to other sets of states, which we briefly summarize here [21] ■ First, the set of classical 
states have the general form 

X= X! p(ki,---,k N )\k 1 ...k N )(k 1 ...k N \. (13) 

ki .. . . ,&n 

This simply means that for a bipartite state of subsystem A spanned by the orthogonal basis {|aj)}, and a subsystem 
B spanned by the orthogonal basis the state x = P{ a i>^j)\ a i)( a i\ ® l&iX&jl is classically correlated when 

one or more of the associated joint probabilities do not factorise as p(a,i, bj) = p{a{) x p(bj). 

The states containing zero correlations, either quantum or classical, are the product states of the form 

7T = 7Tl <g> . . . 7T N . (14) 

Note that the product states are a subset of the classical states which in turn are a subset of the separable states. 
The need to discriminate between separable states and classically correlated states follows from the existence of more 
forms of correlations than just entanglement and classical correlations. This point is discussed further in the following 
section. 



VI. QUANTIFYING DISCORD 

Quantum discord describes all forms of correlations beyond those encapsulated by classically correlated states. For 
a system of two qubits, A and B, the amount discord they share can be quantified using an entropic measure by 
an extension of Shannon's original logic: The more correlated A and B are, the more we can learn about one by 
measuring the other. Accordingly, the degree of correlation between A and B is quantified by the reduction of entropy 
of A (B) when B (A) is measured. 

Measurement of the state of subsystem A is described by a positive operator value measure (POVM) whose elements 
Ej satisfy the completeness relation y\ Ej = 1. For each measurement outcome j, occurring with probability 

Pj = tr (EjPab), the post-measurement state of B is p^ — tiA (E^pab) jVj- The total classical correlations in the 
state pab are then simply the maximum over all possible measurements performed on A of the quantity [2SJ [26] 

C(pab) = S N ( PB ) - 52pjS n (j$). (15) 

3 

This quantity can be considered as one possible quantum generalisation of the classical mutual information Eq. [3] 
The classical information can also be defined by swapping the roles of A and B, but the subtleties related to the 
question of symmetry of correlations in this context will not be relevant for our present discussion. 

The quantum discord of pab is then defined as the difference between the total correlation in the state, given by 
the von Nemann mutual information Eq. ^ , and the classical correlations [23 [21] 

V(B\A) = I(A : B) - max [C\p AB )] 

i E i} 

= min [I (A : B) — C(pab)] , (16) 
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where, as previously mentioned, the maximization is performed over all possible POVMs acting on A. 

The discord is thus the difference between two forms of the mutual information that are equivalent in classical infor- 
mation theory. Hence for classically correlated states the discord is zero, as expected. The discrepancy between these 
two quantities in quantum information theory arises from the role of measurement. While, classically, measurement 
serves to reveal the objective property of the system, in quantum mechanics the outcomes of measurement are basis 
dependent and affect the state of the system. 

More generally, following the discussion of Sec. |Vj for a general N-partite state the discord is quantified by its 
distance to the closest classically correlated state Eq. (13 1. A theory of discord in continuous variable Gaussian states 
has also been developed [21] • Importantly, both entangled and seperable states can have finite discord. Consequently, 
the preparation of states with finite discord is possible under LOCC. This relative ease of preparation has lead to 
questions regarding the technological implications of more general forms of quantum correlations than just entangle- 
ment. Beyond this, the recent level of interest in discord has taught us to discriminate between different forms of 
correlations, a distinction that has no counterpart in classical information. Discord has also received much attention 
in attempts to quantify the "quantumness" of a physical system or process. In the next section we give one such 
application in the context of quantum non-equilibrium thermodynamics. 



VII. CORRELATIONS IN QUANTUM THERMODYNAMICS 

We consider the generation of correlations in two coupled thermal harmonic oscillators subject to a sudden change 
of their harmonic coupling. Starting with a classical Hamiltonian description, we derive the average irreversible work 
done during the change by treating the time dependence of the coupling as a thermodynamic transformation. Next, the 
quantum irreversible work is calculated from an analogous quantised Hamiltonian and transformation. In both cases, 
we assume there is initially zero harmonic coupling and each oscillator is prepared in a Gibbs state, there hence being 
zero correlations between the oscillators. The transformation of the classical system leads to the generation of classical 
correlations only, while, in the quantum case, both classical and non-classical correlations are established between the 
oscillators. We investigate if these more general forms of correlations have a straightforward thermodynamic meaning 
by investigating how the Gaussian discord scales with the difference between the quantum and classical irreversible 
work. For a study of the entanglement properties of coupled thermal oscillators see Refs.|30[ I3T]. 

To begin, we consider two classical harmonic oscillators of identical mass m and natural frequency u) with time- 
dependent coupling X(t), 

rw,\ 1/9 •2^ mUJ 2 . 2 o« wA 2 ' (t) o 

H(X) = — PI + P2) + -iridl + Ql) + -^(li ~ <72) 2 , (17) 
2m 2 2 

where q and p denote position and momentum respectively. Initially, the system is allowed to thermalise with a heat 
bath at temperature /3 = l/fc^T, where fcs is Boltzmann's constant. The partition function of the system is defined 
as E21, 



Z cW ■■= ^ J dqdpe-^W 



2tt\ 2 1 



pojj v /l + 2A 2 > 2 ' 

where q = (gi , 92), P = (jpi , P2) an d h has dimensions of Planck's constant, making the partition function dimensionless. 

The system is prepared by turning off the harmonic coupling (A = 0) and allowing the oscillators to re-thermalise 
before the contact with the heat bath is removed. At t = the system is therefore described by the Gibbs distribution, 

e -0H(o) 

p(q ' p) = ^(op (18) 

where Zc(0) = (27r//i/3wo) 2 is the partition function in the presence of zero harmonic coupling. 

Next, the harmonic coupling is instantaneously switched to the new value A = Xq. In the presence of zero heat flow, 
the first law of thermodynamics asserts that the work done on the system under this protocol is simply the change of 
internal energy 

W c = ff(Ao) - H(0) 
= -~2%i -12) ■ 
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The average work, over many realisations of this protocol, is obtained by averaging Wc over the initial Gibbs distri- 
bution Eq. (fl8|, thus 



(W C ) = J dqdp p(q,p)W, 



c 



dqdp e -/K^M+^^>) 1^° (gi - q2 f 



z c (o) 

\ 2 
f3w 2 ' 

With a little extra effort we are able to calculate the average irreversible (or 'dissipated') work, 

= (W c ) - AF C 
where the free energy change of the process is defined as 

1 Zc(Xo) 



(19) 
(20) 



AF C = --^ln 



P Z c (0) 



2/3 



In 



1 



2A 



2 1 



completing the analysis of the classical system. 

The treatment of the analogous quantum system begins with canonical quantisation of the classical Hamiltonian 
Eq. Jl7l thus 



2m 2 2 



(21) 



where, now, the position and momentum operators obey the canonical commutation relation [ij,pk] = Sj^ih, and ft 
is the reduced Planck's constant. The quantum Hamiltonian Eq. (21) is diagonalised by application of the unitary 
operator (see Ref. [33J for details) 



T — exp 



4ft 



(qtfh -Pih) 



(22) 



Imagining for a moment that X\ t 2 and pi 2 represent the quadratures of two modes of the electromagnetic field, the 
operator T has the physical interpretation of mixing the two modes at a 50:50 beam-splitter. Under this operation, 
the Hamiltonian Eq. (21) becomes 



H' = THT^ 



1 



Pi 



muj\ 2 

2m- 1 ' ~ Ql + 2^ Pl ' 2 

where the harmonic coupling is encoded in the renormalized frequencies wi(A) = WoyT — A 2 /wq, UJ2(X) 
ujq^I + X 2 /ujq and uj = \/u> 2 + A 2 , thus wi(0) = u; 2 (0) = ui. 
Finally, defining the creation and annihilation operators 



1 . 2 tow 2 2 



"92, 



(23) 



flfc(A) := 



4(A) := 



2ft 
2ft 



<7fc H Pfe 

TOWfe 



TOWfe 



that obey the commutation relation 



= 5j k VA, Eq. (231 can be rewritten as 



H' = hui ( a\ai + -z) + f^2 [ a\a 2 + 



1 



(24) 



(25) 



Eqs. (23) and (25) have the familiar form of two uncoupled harmonic oscillators with natural frequencies uj\ and u)2 

/erse t( 

-pit' 1 , fPtuJi 



respectively. Consequently the partition function of the system at inverse temperature j3 is simply the product of the 
partition functions for each individual oscillator i.e., 



Z Q (X) :=tr 



-PH 



= tr 



= -csch 
4 



csch 



where tr[.] denotes the trace over the two-mode Fock basis spanning the Hilbert space of the Hamiltonian. 

To proceed, we consider an identical protocol to that performed upon the classical system: Initially, the harmonic 
coupling is zero (A = 0) and the system is prepared in the corresponding Gibbs state at inverse temperature /3, 



1 



PG 



Zq(0) 



-/9H(0) 



(26) 



where the appropriate partition function is explicitly Zq(0) = (l/4)csch 2 (/3/ia;/2). Subjecting the oscillators to an 
identical thermodynamic transformation, taking the harmonic coupling to the final value A 
following object, 



An, we define the 



W Q = H(Xa) - H(0) 



<12 ) 



HX 



°(a t 1 (0) + a t 1 (0)-a 2 (0)-aj(0)) 2 



(27) 



Note that, in general thermodynamic transformations of quantum systems, it is not possible to define a so-called 
'operator of work' |34j . However for a sudden quench of the harmonic coupling, Eq. (27) allows definition of the 
average work over many realisations of the protocol as, 



{Wq) = tr [ Pg Wq 

i h\l 



Zq(0) 2uj 

nxl , 

= TT^ COth 

2w 



tr 



-^(al(0)a 1 (0)+ala 2 +i) (al (0)ax (0) +a f 2 a 2 



(28) 



Again, we are able to calculate the average irreversible work, 



(Wq") = (Wq) - AFq, 



(29) 



where the free energy change is, 



1 ln Zq(X ) 
P Zq(0) 

1 l csch(/3fiwi(A )/2)csch(/3fio; 2 (Ao)/2) 
n csch 2 (jSfiw/2) 



(30) 



The thermodynamic significance of the the quantities Eqs. (28) - &30h can be understood as follows: For a quasistatic 



transformation, the irreversible work is zero and the total work is given simply by the free energy change. The free 
energy change is thus the energy required to transform between spectra of the initial and final Hamiltonian. For 
non-quasistatic processes, however, the average irreversible work is finite and gives the energy expended in exciting 
the system during the transformation. 



In the quantum regime ksT << ftwn the expressions for the quantum and classical irreversible work (Eqs. (29 1 and 



(20), respectively) differ. Wc thus define the average excess quantum dissipated work as 



n 



Q 
ftA 2 
~2uj 



diss^ 



coth 



/3huj 



(3huj 



1 csch ((3huj 1 (X )/2) csch (f3tuv 2 (X )/2) 
P n csch 2 {phu/2) 



2(3 



In 



2A 2 



Clearly, fl > 0, with equality reached in the classical limit ksT » hujo or for a quasistatic process. The average 
excess quantum dissipated work is thus the additional energy expended in exciting a quantum system relative to its 
classical analogue undergoing an identical transformation. 

We investigate whether this quantity also encodes the fact that during the transformation, the two components of 
the quantum system share additional, more general forms of correlations than available to its classical analogue. As a 
starting point, following Ref. |29j we calculate the Gaussian discord induced by the sudden switch on of the harmonic 
coupling. Explicitly, we consider the discord T> generated in transforming an initial product of identical thermal states 
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FIG. 1: The average excess quantum dissipated work (dotted line) and quantum discord (solid line) following a sudden quench 
of the harmonic coupling with amplitude Ao = ui and h = kB=u = m = l, as a function of temperature. 



Eq. (26) according to the unitary evolution U — e jH ( A °' generated by the Hamiltonian Eq. (21 1 i.e., T>(U pcU^). Note 
that as the Hamiltonian is bilinear in the creation and annihilation operators, it preserves the Gaussian character of 
the initial thermal states. 

At this stage, without loss of generality we set ks = h = u> = m = 1, and plot (Fig. [T]) the average excess quantum 
work and Gaussian discord for a quench amplitude Ao = u as a function of temperature T. Fig.[l]shows that the discord 
and excess quantum irrevesible work show qualitatively similar behaviour, decaying with increasing temperature as 
expected. However, it is not clear if there exists a definitive link between the two quantities. While the excess quantum 
irreversible work should contain some information on the thermodynamic cost of correlation generation, this is likely 
to be convoluted with contributions to the total from other sources, e.g. exciting the oscillators individually without 
correlating them. Nonetheless, it remains an interesting open question as to what the relation of quantum discord, 
and other forms of correlation, to non-equilibrium thermodynamics may be. 
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